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NOTATION 


cylindrical  coordinates 

inner  and  outer  radii,  respectively,  of 
body  of  revolution 

ply  thickness 

involute  constant 

number  of  plies 

angle  of  rotation  of  EIS 

arc  angle 

tilt  angle 

surface  angle 

helical  angle 

meridian  angle  in  ply  pattern 

helical  convolute  angle 

constant  defined  by  Eq.  (40) 

cylindrical  involute  constant 

helical  arc  length  of  helical  convolute 

tangent  length  of  helical  convolute 

radius  of  cylinder  in  ply  pattern 

arc  angle  in  ply  pattern 

plane  polar  coordinates  in  ply  pattern 

minimum  and  maximum  values  of  Z  in  EIS 

distance  between  reference  points  in 
ply  pattern 

surface  area  of  EIS 

arc  length  along  involute 

volume  of  body  of  revolution 


in  > 


1,  ],  k  =  unit  vectors  along  0,  r,  z  directions, 

respectively 

,  n  =  unit  vectors  on  EIS 


Overbars  stand  for  parameters  at  arbitrary  fixed  point  P  on  EIS. 
Primes  stand  for  undebulked  parameters  or  0'  meridian. 


INTRODUCTION 


Involute  construction,  formerly  known  as  rosette  construction, 
is  a  popular  approach  being  used  in  rocket  nozzle  technology  for 
the  fabrication  of  exit  cones  and  other  bodies  of  revolution. 

Such  bodies  are  formed  by  laminating  identical  fabric-reinforced 
composite  plies  of  uniform  thickness  in  such  a  way  that  each 
ply  extends  to  the  extremities  of  the  body  in  both  the  radial 
and  axial  directions.  Hence,  the  boundaries  of  the  body  are 
completely  defined  by  the  ply  edges.  The  mathematical  solution 
of  this  problem  is  the  theme  of  this  work. 

The  contemporary  approach  to  involute  construction  incorpor¬ 
ates  the  use  of  analysis  and  fabrication  techniques  that  can 
lead  to  significant  pre-existing  flaw  distributions  within  the 
body.  Defects  such  as  wrinkles  and  fabric  distortion  are  induced 
and  represent  the  suspected  source  of  numerous  failures  that 
have  occurred  in  test  firings,  and  even  in  fabricated  articles 
prior  to  imposed  loadings.  Regions  of  geometrical  discontinuity, 
such  as  the  neighborhood  of  a  cylinder-cone  intersection  are 
particularly  sensitive  to  the  formation  of  this  initial  damage. 

The  origin  of  many  of  these  problems  appears  to  lie  in  the  use 
of  an  incorrect  ply  pattern,  i.e.,  the  basic  ply  geometry 
adopted  does  not  satisfy  the  involute  construction  problem,  which 
requires  that  a  prescribed  volume  of  revolution  be  filled  pre¬ 
cisely  with  a  given  number  of  identical  plies  of  the  same  constant 


thickness. 
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The  first  analytical  treatment  of  involute  bodies  appears 
to  be  due  to  Mamrol  (Ref.  1)  in  an  unpublished  report.  Other 
unpublished  work  was  accomplished  by  workers  at  Hitco  Corp. 

(Refs.  2,3).  Pagano  (Ref.  4)  provided  the  first  comprehensive 
treatment  of  cylindrical  involute  bodies,  including  the  appro¬ 
priate  governing  equations  for  the  stress  field  in  an  elastic 
body  under  axisymmetric  loading.  Subsequently,  an  approximate 
method  was  presented  (Ref.  5)  to  establish  the  ply  pattern  for 
bodies  of  conical  shape.  In  the  latter  reference,  it  was  assumed 
that  one  edge  of  the  ply  pattern  is  a  straight  line  that  coincides 
with  the  generator  of  a  cone.  It  will  be  shown  here  that  this 
assumption  is  only  valid  under  very  restrictive  conditions.  In 
most  cases,  the  work  in  the  study  of  involute  bodies  proceeds 
by  first  establishing  the  theoretical  surface  (called  the 
involute  surface)  that  satisfies  the  involute  construction 
problem,  following  which  an  approach  to  develop  that  surface, 
i.e.,  the  ply  pattern  generation  scheme,  is  presented.  The 
inverse  problem  of  definition  of  the  involute  surface  given  a 
specific  ply  pattern  was  treated  by  Stanton  and  Pagano  (Ref.  6) , 
where  an  axisymmetric  stress  field  model  was  also  derived  and 
used  to  examine  the  tendency  of  certain  exit  cones  to  delaminate 
under  processing  conditions.  Since  that  treatment  of  the 
involute  surface  geometry  was  also  approximate,  the  question  of 
existence  of  a  solution  to  the  involute  construction  problem  was 
naturally  raised.  In  this  work,  we  shall  address  this  question. 


2 


In  a  parallel  effort  with  the  present  work,  E.E.  Savage 
(Ref.  7)  has  postulated  that  a  helical  convolute  surface  satisfies 
the  requirements  of  the  involute  construction  problem.  We  shall 
provide  the  proof  of  this  assertion  in  this  presentation,  as  well 
as  to  rigorously  derive  the  complete  mathematical  formulation 
that  characterizes  the  aforementioned  solution,  including  the 
equations  that  describe  the  ply  pattern.  An  experimental  demon¬ 
stration  of  the  present  formulation  shall  also  be  given,  as  well 
as  a  model  to  characterize  the  kinematical  response  during  the 
debulking  process,  which  is  an  important  part  of  the  fabrication 
method . 


DEFINITION  OF  EXACT  INVOLUTE  SURFACE  (EIS) 


Consider  a  body  of  revolution  generated  by  rotation  of  a 
plane  region  about  the  z-axis  as  shown  in  Figure  1.  The  boundary 
of  the  plane  region  is  denoted  by  C,  and  conventional  cylindrical 
coordinates  (r,  9,  z)  will  be  employed  in  the  formulation.  Left- 
and  right-hand  segments  of  the  contour  C  are  introduced  and 
defined  by  radial  coordinates  r  =  R^  (z),  (z ) ,  respectively. 

These  two  segments  may  either  constitute  the  entire  boundary  C, 
as  in  Figure  1(a),  or  be  separated  at  the  upper  and/or  lower 
extremities  by  straight  lines  normal  to  the  z-axis,  such  as  the 
case  illustrated  in  Figure  1(b).  In  general,  the  curve  C  must 
be  represented  by  a  piecewise-continuous  function. 

In  practice,  involute  const ruction  of  bodies  of  revolution 
such  as  those  shown  in  Figure  1*  employs  a  lamination  procedure 
utilizing  identical  plies  of  a  constant  thickness.  This  practice 
leads  to  the  definition  of  the  involute  construction  problem, 
which  can  be  posed  as  a  hypothetical  problem  that  admits  to  a 
mathematically  exact  treatment,  i.e.,  determine  the  equations  of 
identical  surfaces  which  are  separated  everywhere  by  a  constant 
distance  T  measured  in  the  local  normal  direction  to  the  surface 
such  that  their  boundaries  completely  define  the  surface  of  a 
given  body  of  revolution  (see  Figure  2) .  j 

Hence,  the  volume  defined  by  the  boundaries  of  surfaces  S  j 

j 

is  equivalent  to  the  volume  inside  the  body  of  revolution.  If  j 

1 

the  body  were  actually  formed  from  material  layers  or  plies  of  i 

J 

thickness  T  whose  boundaries  follow  the  hypothetical  surfaces  S,  i-j 


Figures  are  located  in  Appendix  C 
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interior  of  the  body  would  be  filled  such  that  no  gaps  occur 
between  contiguous  plies.  A  slight  difference  would  exist  between 
the  material  volume  and  the  theoretical  volume  of  the  body  of 
revolution  because  of  the  step- like  appearance  created  by  the 
plies  at  the  boundaries.  However,  this  difference  vanishes  as 
T-K>  and  does  not  exist  at  all  in  the  hypothetical  problem.  The 
solution  of  the  mathematical  analog  of  the  physical  problem  will 
be  referred  to  here  as  the  exact  solution  and  the  surface  S  will 
be  called  the  exact  involute  surface  (EIS) . 

TRANSFORMATION  RELATIONS 

In  this  section,  we  shall  formulate  the  approach  which 
leads  to  the  equation  of  the  EIS.  Our  approach  consists  of 
examination  of  the  relationship  that  exists  between  two  surfaces, 

S  and  S',  within  the  given  body  of  revolution.  The  two  surfaces 
are  identical  in  size,  shape,  and  all  local ■ geometric  details, 
except  that  they  are  displaced  from  one  another  in  the  body. 

We  shall  discuss  two  successive  transformations,  which  carry  S 
into  S ' ,  thence  S '  into  S ,  such  that  the  requirements  of  the 
mathematical  analog  are  satisfied  (see  Figure  3) . 

In  the  first  transformation,  each  point  in  S  translates 
through  the  displacement  Tn,  where  n  is  the  local  unit  normal 
vector.  This  transformation  maps  S  into  S'. 

The  second  transformation  follows  from  the  consideration  of 
the  geometry  in  the  cylindrical  involute  problem  (Ref.  4) .  In 
that  case,  we  have  shown  that  an  annular  region  can  be  filled 
with  identical  curves,  called  involutes  or  spirals,  which  are 
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related  to  one  another  by  rigid  body  rotations  about  the  z-axis. 
Since  each  cross-section  of  our  body  of  revolution  is  an  annulus, 
this  transformation  is  also  appropriate  in  the  present  example. 
Hence,  our  second  transformation,  which  carries  S'  back  into  S, 
consists  of  a  rotation  through  an  angle  $  about  the  z-axis,  where 
$  is  a  constant. 

The  two  transformations  above  guarantee  that  the  volume  of 
a  body  of  revolution  can  be  filled  with  identical,  uniformly 
spaced  surfaces.  As  the  surface  S  may  extend  indefinitely,  we 
must  truncate  it  at  its  intersection  with  the  surface  produced 
by  revolution  of  curve  C. 

MATHEMATICAL  FORMULATION 

Consider  the  trace  of  surface  S  that  lies  in  the  plane  z  = 
constant.  It  has  been  shown  (Ref.  4)  that  a  plane  annulus  can 
be  filled  by  identical  OUrves  (involutes)  related  to  each  other 
by  rigid  body  rotations  about  z.  These  curves  are  defined  by 
the  equation 

r  sin  a  =  c  (1) 

where  c  is  a  constant  and  a,  called  the  arc  angle,  is  the  angle 
between  tangents  to  the  curve  and  the  circle  of  radius  r  at  the 
same  point  on  the  curve.  In  order  to  satisfy  our  second  trans¬ 
formation,  an  equation  of  the  form  (1)  must  be  satisfied  within 
each  cross-section  of  a  general  body  of  revolution.  Thus,  the 
general  solution  for  the  EIS  takes  the  form 

r  sin  a  =  c  (z)  (  2) 

where  c  is  a  function  of  z  alone.  It  will  now  be  necessary  to 
review  some  of  the  basic  relations  needed  to  characterize  sur¬ 
face  geometry  (Ref.  6) . 


•  r 
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Consider  an  infinitesimal  element  of  a  surface  as  shown  in 


Figure  4 .  Two  right-handed  systems  of  unit  vectors  are  shown 
in  the  figure:  first,  letting  carets  symbolize  vector  quantities, 
we  have  the  orthogonal  system  i,  3»  £»  which  lie  in  the  positive 
0,  r,  and  z  coordinate  directions,  respectively;  secondly,  the 
non-orthogonal  system  t,  s,  n.  Here  t  is  the  vector  tangent  to 

A 

the  meridian  0  =  constant,  s  is  tangent  to  the  trace  of  the 

A 

surface  in  the  cross-section  z  =  constant,  and  n  is  the  normal 
vector  pointing  away  from  the  interior  of  the  body.  Arc  angle  a, 
tilt  angle  y,  and  surface  angle  ip  are  also  shown  in  Figure  4,  where 

tan  Y  =  af  (3) 

A  A 

cos  ip  =  t  •  s  (4) 


We  are  now  ready  to  invoke  the  first  transformation  of  the 
preceeding  section,  i.e.,  local  translation  along  n. 

We  start  from  an  arbitrary  point  p  (r,  6,  z)  that  lies  in 
surface  S.  From  this  point,  we  construct  a  normal  vector  ATn, 
which  intersects  S'  at  p'  (r  +  Ar,  6  +  A6,  z  +  Az)  and  is  normal 
to  S'  at  the  latter  point.  Letting  ft'  represent  the  position 
vector  of  p ' ,  we  see  that 

R'  =  r3  +  ATn  (11) 

Substitution  of  Eq.  10  into  (11)  leads  to 

(r  +  Ar)  2  sin2  ip  =  r2  sin2  ip  +  (AT)2  cos2  y  (12) 

+  2rAT  sin  cos  a  cos  y 


along  with 


and 


sin  A9  =  sing  cos  y 
(r  +  Ar)  sm  ^ 


Az  =  -AT  cos  a  sin  y 
sin  if 


(13) 


(14) 


We  may  now  write  an  expression  for  n',  the  unit  outward  normal 
vector  at  p1  in  terms  of  a',  y',  and  '  ,  and  the  basis  vectors 
£ ' ,  j ' ,  and  k ' ,  where 

i'  =  i  cos  A0  -  3  sin  A0  (15) 

j"  =  i  sin  A0  +  j  cos  A0  (16) 

Substituting  Eqs.  (15)  and  (16)  into  this  expression  and  using 
the  fact  that  n'  =  n,  we  can  show  that 

g'  “  g  =  Ag  =  “A0  *  (17) 

where  a'  is  the  arc  angle  at  p'. 
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Let  us  next  express  Eq.  (2)  at  point  p',  with 


Ac  =  c  (z  +  az)  -  c  ( z )  (18) 

so  that 

(r  +  Ar)  sin  (a  +  Aa)  =  c  +  Ac  (19) 

Inserting  Eq.  (17) ,  this  becomes 

(r  +  Ar) (sina  cosA9  -  cosasinA9)  =  c  +  Ac  (20) 

However,  because  of  Eq.  (13) ,  we  find  that 


_  r  sm  +  AT  cos  y  cos  a 

COS  At)  =  - 5 - ; - 1 - r - r-i - ; - 

(r  +  Ar)  sin  ip 


(21) 


Finally,  putting  Eqs.  (13)  and  (21)  into  (20)  yields 

r  sin  a  =  c  +  Ac  (22) 

which,  along  with  Eq.  (2) ,  implies  that  Ac  =  0,  or 

r  sin  a  =  c  =  constant  (23) 

Hence,  the  EIS  must  satisfy  Eq.  (23)  in  order  to  obey  the  first 
required  transformation  (ATn) . 

We  shall  now  seek  to  establish  the  appropriate  relations 
such  that  the  second  transformation  (rigid  body  rotation)  is 
satisfied.  We  first  observe  that  Eq.  (23)  implies  the  following 
relation  for  a  given  value  of  z , 

9(a,z)  -  ¥(2)  =  a(z)  +  cot  a(z)  -  a(z)  -  cot  a(z)  (24) 

since  Eq.  (23)  is  the  same  equation  that  governs  cylindrical 
involute  geometry  (Ref.  4).  Here,  P(r(z),  ^(z) ,  z)  represents 
an  arbitrary  fixed  point  in  S,  while  ct(z)  is  the  arc  angle  at 
this  point,  and  we  have  chosen  to  regard  9  as  a  function  of  a  and 
z  in  this  relation.  Of  course,  a  and  r  are  related  everywhere  via 
Eq.  (23) .  For  conceptual  reas  ns,  one  may  choose  the  arbitrary 
fixed  point  to  coincide  with  the  intersection  of  the  EIS  at  z 
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with  the  inner  or  outer  surface  of  the  body  of  revolution,  i.e., 
R^(z)  or  I^f2)  as  s^own  Figure  1,  although  this  is  not 
required. 

In  order  to  treat  the  rigid  body  rotation,  it  is  necessary 
to  consider  the  traces  of  two  involute  surfaces,  S  and  S',  cut 
by  a  plane  z  =  constant  since  these  traces  map  into  each  other 
in  such  a  transformation.  As  before,  we  begin  by  constructing 
a  vector  ATn  from  point  P  such  that  the  tip  of  the  vector  inter¬ 
sects  S'  at  P'  (r(z)  +  Ar,  7T(Z)  +  A"0,  z  +  AZ)  .  Using  primes  to 
denote  functions  in  S'  and  taking  Eq.  (17)  into  account, 

Eq.  (24)  gives 

9'  (a',  z  +  az)  -  (0(z)  +  a'S)  =  ct(z)  -  Al  +  cot(a(z)  -  A9) 

(25) 

-  a' (z  +  Az)  -  cot  a' (z  +  Az) 
Equation  (24)  with  z  replaced  by  z  +  az  and  Eq.  (25)  are  the 
equations  of  the  traces  of  S  and  S',  respectively,  at  the  level 
z  +  AZ.  The  rigid  body  rotation  is  therefore  expressed  by 

0 '  (a,  z  +  AZ)  -  0(a,  z  +  Az)  =  A<D  (26) 

since  points  with  the  same  a  also  have  the  same  r.  In  Eq.  (24) , 

A4>  is  a  constant  to  be  evaluated  later  on.  Carrying  out  Eq.  (26)  , 
with  the  aid  of  (24)  and  (25),  leads  to 

A$  +  AIT  =  a(z)  +  cot  a(z)  -  a(z  +  Az)  -  cot  a(z  +  AZ) 

Az  cot  y(z) 
c  cos  a(z) 

where  Eqs.  (13),  (14),  and  (21)  have  been  used,  as  well  as 
elementary  trigonometric  formulas.  Also,  we  have  defined 

&W  =  e^z  +  Az)  -  lf(z) 
and  y(z)  is  the  tilt  angle  y  at  P. 


(27) 


(28) 


In  order  to  solve  Eq.  (27)  for  AF,  we  only  need  to 
determine  the  constants  A4>  and  y(z)  in  terms  of  known  parameters 
since  a  (*)  is  known  from  Eq.  (23)  if  r(z)  is  given,  e.g.,  R^(z) 
or  R2 (a) . 

Considering  the  angle  A4>  first,  we  note  that 


A$ 


2tt 

AN 


(29) 


where  AN  is  the  number  of  "layers"  or  exact  involute  surfaces 
corresponding  to  thickness  AT  required  to  fill  the  annular  space. 
Hence,  we  get 

NT 

AN  =  (30) 

where  N  is  the  number  of  "layers"  of  thickness  T.  Therefore 


A  <S>  = 


2itAT 

NT 


which,  by  virtue  of  Eq.  (14) ,  becomes 

A$  =  -2ttAz  sin  \jT  (z) 

NT  cosci(z)  siny(z) 


(31) 


(32) 


where  ty(z)  is  the  surface  angle  ip  at  P. 

To  define  the  tilt  angle  y(z),  we  consider  the  meridian  of 

the  involute  surface  S  containing  point  P.  Let  this  meridian 

intersect  the  z  +  Az  trace  of  S  at  r  =  r  (a  -  a  ) .  Thus,  an 

P  P 

approximate  expression  for  y(z)  is  given  by 


r  -  r(p) 
tan  y  (z )  =  — *- 


Az 


(33) 


which  becomes  exact  in  the  limit  as  Az-*0. 

Writing  Eq.  (24)  at  z  +  Az  and  putting  0 (a  ,  z  +  Az)  =  F(z) 
yields  the  result 


A0  =  a  +  cot  ci  -  a(z  +  Az)  -  cot  a(z  +  Az) 
P  P 


(34) 
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(35) 


which  brings  Eq.  (27)  to 

A4>  =  a(z)  +  cot  ct(z)  -  a  -  cot  a  -  -z--cot 

”  P  c  cos  a(z) 

Now  using  Eqs.  (23)  and  (33)  and  elementary  trigonometric 


identities,  we  can  show  that 

tan  y(z)  tan  a(z)  Az 
r  (z) 

in  the  limit  as  Az+0,  so  that 


ap  =  a(z)  - 


,  — ,  .  .  tan  y(z)Az 

cot  a  =  cot  a(z)  +  - — — 


(36) 


(37) 


c  cos  a(z) 


Finally,  substituting  Eqs.  (32),  (36),  and  (37)  into  (35)  while 
passing  to  the  limit  Az+0,  we  find  that 


2iTsinjMz) 


~  tan  y(z)  +  cot  y(z) 


NT  cos  a(z)  sin  y(z)  c  cos  a(z) 


c  cos  a(z) 


(38) 


tan  a(z)  tan  y(z) 
r  (z) 

which,  after  some  elementary  manipulations  and  the  use  of  Eq.  (7) , 


yields 


where 


tan  y(z)  = 


cos  a(z) 


A  =  — ^  (4irzc2  -  N2T2)  ** 
which  is  a  constant.  Clearly,  we  have 

c  > 


(39) 


(40) 


(41) 


which  controls  the  minimum  value  of  c  that  can  be  utilized. 

Given  the  expressions  (32)  and  (39),  we  may  now  return  to 
Eq.  (27)  to  determine  A¥.  We  now  define 

Aa  =  a(z  +  Az)  -  a(z)  (42) 

Introducing  Eqs.  (39),  (42),  and  (32),  in  conjunction  with  (7), 
into  (27)  and  letting  Az->-0,  we  arrive  at 
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n 


de  =  £  dz  +  cot2  a(z)  da  (43) 

c 

On  integration  over  the  range  (zq,  z) ,  Eq.  (43)  gives 

1)(z)  -  0(z  )  =  ^  (z  -  z  )  +  a(z  )  +  cot  a (z  )  -  a(z) 

o  c  o  o  o  (44) 

-  cot  a(z) 

and  by  substitution  into  Eq.  (24) ,  we  get 

0(a/Z)  -  "^(z0)  =  §  <z  “  zq)  +  «(z0>  +  cot  o(zQ)  -  a (z)  ^ 

-  cot  a(z) 

which,  along  with  (23)  ,  represents  the  general  equation  of  the 
EIS.  The  point  corresponding  to  a  =  a(zQ)  ,  0  =  ^(zq)  =  0(a,  zq) , 
z  =  zQ  simply  represents  an  arbitrary  point  on  the  EIS.  Also, 
since  this  point  is  arbitrary,  Eq.  (39)  can  be  written  as 


tan  y  (a,  z)  = 


cos  a(z) 


so  that  the  relation  (46)  is  valid  everywhere  on  the  EIS. 

An  important  feature  inherent  in  the  present  formulation  is 
the  clear  statement  of  the  "input"  or  prescribed  parameters 
necessary  to  define  the  EIS.  Any  combination  of  parameters 
that  serve  to  uniquely  define  the  following  constants  and 
functions  constitute  a  necessary  and  sufficient  input  data  set, 

c,  N,  T,  R^z)  ,  R2  (z)  ,  0Q  (47) 

where  0Q,  the  0-coordinate  of  any  point  in  R1(z),  or  R2(z),  is 
only  needed  to  define  the  angular  position  of  the  EIS,  not  its 
shape  or  size,  and  R1(z)  and  R2(z>  need  only  be  piece-wise 
continuous  functions.  While  contemporary  manufacturing  procedures 
and  the  approximate  treatment  (Ref.  5)  assume  that  a  trace  of 
the  involute  surface  in  a  meridional  plane  coincides  with  the 
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generator  of  a  cone,  the  input  parameters  (47)  ,  which  define 
the  EIS,  prohibit  this  assumption,  in  general.  However,  if  we 
invoke  the  small  a  assumption  a2<<a,  we  find  that  Eq.  (44)  is 
satisfied  by  the  conical  trace  r(z)  =  kz  provided  that  A  =  K, 
where  A  is  given  by  Eq.  (40)  and  K  is  a  constant.  In  this  case 
then,  the  conical  trace  is  a  good  approximation. 

Apparent  singularities  occur  in  Eqs.  (27) ,  (32) ,  (35) ,  and 
(38)  for  the  case  y(z)  =0.  These  are  removed  by  appealing  to 
Eqs.  (7)  and  (14) ,  which  give,  for  example 


-AT  cos  a 


(48) 


The  situation  y(z)  =  0  arises  only  when  the  constant  A  vanishes. 


This  leads  to 


(49) 


and  corresponds  to  a  cylindrical  EIS,  which  is  the  significance 
of  subscript  c  in  Eq.  (49) .  The  present  equations  now  reduce 
identically  to  those  given  earlier  (Ref.  4)  for  cylindrical 
involute  surfaces.  In  this  regard,  it  is  important  to  note  that 
the  generation  of  a  cylindrical  body  does  not  require  that  the 
involute  surface  be  cylindrical.  An  example  of  this  will  be 
shown  later  here. 


The  surface  given  by  Eqs.  (45)  and  (23)  represent  the  exact 
solution  for  the  mathematical  problem  of  filling  space  with 
identical,  uniformly  spaced  surfaces.  These  equations  represent 
a  helical  convolute  surface  (Ref.  8) .  In  the  next  section, 


we  shall  examine  some  of  the  basic  parameters  of  the  helical 
convolute  since  they  lead  to  an  expedient  approach  to  develop 
the  surface,  i.e.,  to  construct  the  ply  pattern. 


DEVELOPMENT  OF  EIS 


We  now  proceed  to  demonstrate  that  the  EIS  is  a  helical 
convolute  of  a  cylinder  of  radius  c.  To  start,  we  shall  review 
the  geometrical  construction  of  the  helical  convolute  (Ref.  8) , 
which  will  be  helpful  in  establishing  the  equation  of  the 
developed  surface,  or  ply  pattern. 

The  helical  convolute  is  generated  by  a  straight  line  which 
moves  in  space  such  that  it  remains  tangent  to  a  helix.  The  case 
where  the  end  point  of  the  line  moves  in  a  plane  normal  to  the 
axis  of  the  cylinder,  which  is  important  in  the  construction  of 
involute  exit  cones;  is  depicted  in  Figure  5.  In  this  figure,  the 
helix  on  the  cylinder  of  radius  c  is  shown  as  OQT.  The  helical 
convolute  is  generated  by  PQO  such  that  the  sum  of  the  helical  arc 
length  w  =  OQ  and  the  tangent  length  l  =  QP  remains  constant. 

In  other  words,  the  helical  convolute  is  generated  by  wrapping  a 
right  triangular  sheet  of  height  OS  and  hypotenuse  PQO  on  the 
cylinder.  As  a  result  of  this  geometrical  construction,  we  can 
show  that  the  coordinates  of  points  P  and  Q  are  related  by 


r2  =  (z  -  Zg) 2  tanz6  +  c2 


sin  (n/2  -  0  +  0q)  =  c/r 


(50) 


By  direct  substitution,  one  can  demonstrate  that  the  relations 
(50)  identically  satisfy  the  EIS  equations  (45)  and  (23) , 
provided  that 

tan  6  =  A  (51) 

Now  the  angle  5  is  defined  by 


tan  6  = 


cd0 


dz 


(52) 
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If  we  let  the  helix  OQT  represent  the  intersection  of  the  EIS 

and  the  cylinder  of  radius  c,  Eq.  (45)  gives 

d0  =_A_dz  (53) 

c 

Substituting  Eq.  (53)  into  (52)  shows  that  Eq.  (51)  is  satisfied. 
Hence,  the  postulate  that  the  EIS  takes  the  form  of  a  helical 
convolute  (Ref.  7)  has  been  proven.  Finally,  the  tangent  length 
i  is  given  by 


l  -  (z  -  z q)  /  1  +  A2  (54) 

In  order  to  employ  the  more  convenient  base  point  0,  it  is  only 
necessary  to  substitute  the  coordinates  of  point  0  in  place  of 
those  for  Q  in  the  above  relations. 

Since  the  EIS  is  a  helical  convolute,  it  represents  a 
single-curved,  developable  surface  (Ref.  8) .  It  follows  that  the 
helix  of  radius  c  is  mapped  into  the  arc  of  a  circle  of  radius 
Rq,  where 


Ro  sin*  5 


(55) 


which,  because  of  Eq.  (51) ,  becomes 

Rq  =  c  (1  +  —>)  (56) 

Also,  the  lengths  w  and  l  (Figure  5)  are  preserved  in  the  mapping 
so  that  the  path  of  any  point  such  as  P  (which  is  an  involute 
of  a  circle  of  radius  c)  is  carried  into  an  involute  of  a  circle 
of  radius  RQ.  Reference  to  Figure  6  indicates  that  the  mapped 
image  of  point  P  is  described  by  polar  coordinates  R,  ®  ,  where 

®  ■  tt/2  +  cot  XQ  -  X  -  cot  X  (57) 


R  sin  X  =  R, 


(58) 
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and 


tan  X 


(z  -  z„J  /  1  +  A2 


Using  the  first  of  Eqs.  (50) ,  along  with  (51)  and  the  relation 
that  derives  from  the  preservation  of  arc  length  w,  namely 


ZQ  '  Zp 


=  Ro  (  @  +  X  -  tt/2) 
/  1  +  A2 


we  find  that  Eq.  (57)  can  be  put  into  the  form 


cot  X= 


/TTT2  ( 

Ro  V 


-  z  +  z  1  +  cot  X 


and  finally,  using  Eqs.  (59)  and  (56),  we  get 


cot  X  =  — 

c 


' r2  -  c2 
1  +  A2 


which  can  also  be  expressed  as 


cot  X  = 


/l  +  A2 


cot  a 


Hence,  once  the  EIS  has  been  defined,  an  arbitrary  point  (r,  0, 
z)  in  space  is  mapped  into  (R,  (g)  )  in  developed  space  by  solving 
Eq.  (62)  for  X  and  substituting  into  (57)  and  (58) ,  after  using 
(56)  and  (59)  to  compute  XQ  and  Rq.  It  must  be  emphasized  that 

we  have  implicitly  assumed  that  z  -  z  >  o  so  that  z  must  be 

p  —  p 

chosen  as  the  minimum  value  of  z.  The  value  of  z^  corresponds  to 
the  point  0  in  Figure  5  at  which  w  =  0  and  is  easily  found  from 

the  equations  of  the  EIS  by  repeated  use  of  the  relation 

y? z  _  _2 

zp  =  Z  -  *4-*  (64) 


Here  r  and  z  are  the  coordinates  of  any  point  on  the  EIS  and 


Eq.  (64)  is  used  to  search  for  the  minimum  value  of  z^,  which 
will  always  correspond  to  a  point  on  the  outer  boundary  of  the 
body  of  revolution  being  generated.  Finally,  since  the  EIS  is 
a  developable  surface,  we  shall  dispense  with  the  distinction 
between  the  mathematical  and  physical  solutions  of  the  involute 
construction  problem;  they  only  differ  because  of  the  finite 
compliance  of  the  sheet  material  and  the  "steps"  on  the  boundary 
of  the  body  constructed  of  sheets  of  constant  thickness,  as 
discussed  earlier. 


ILLUSTRATIVE  EXAMPLE 

As  an  example  of  the  use  of  the  present  model,  the  body 
of  revolution  Having  the  profile  shown  in  Figure  7  was  built  by 
involute  construction.  The  geometry  of  this  body  is  similar  to 
that  of  contemporary  rocket  nozzle  exit  cone  structures.  The 
ply  material  used  in  the  construction  was  heavy  paper  having  a 
constant  thickness  of  0.01  in.  Using  this  value  for  T,  with 
c  =  0.15  and  N  =  92,  which  satisfy  the  inequality  of  Eq.  (^.1)  , 
the  ply  pattern  generated  via  the  present  model,  Eqs.  (56)— (59) , 
(62),  and  (64),  is  shown  in  Figure  8.  The  mapped  positions  of 
several  involute  trajectories  (z  =  constant)  and  meridians 
(0  =  constant) ,  which  approximate  straight  lines  for  reasonably 
small  values  of  a,  are  also  shown  in  the  figure.  The  range  of 
Y  in  this  body  is  given  by  12.56°  <_  y  <  12.78°. 

Following  the  ply  pattern  generation,  the  body  of  revolution 
was  constructed  by  laminating  92  plies  such  that  their  inner 
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boundaries  rested  on  the  surface  of  the  male  mandrel  defined  by 
radii  r  =  R^(z).  Pre-forming  the  plies  is  convenient,  though 
not  necessary.  A  pre-formed  ply  in  the  shape  of  the  EIS  and 
the  male  mandrel  are  shown  in  Figure  9,  as  well  as  the  body 
itself.  In  practice,  the  plies  would  be  bonded  together,  however, 
in  the  present  case,  they  are  held  in  place  by  supporting  rings 
as  shown  in  Figure  9.  The  existence  of  this  body,  in  accordance 
with  the  geometry  prescribed  in  Figure  7 ,  serves  as  an  experimental 
demonstration  of  the  validity  of  the  present  mathematical  model. 

The  body  shown  in  Figure  9  could  not  be  built  using  contemporary 
practice  since  the  heavy  paper  material  has  no  mechanism  to 
produce  the  severe  distortions  associated  with  that  practice. 

Thus  the  potential  to  build  relatively  defect-free  structure  is 
another  virtue  of  the  new  approach. 

There  exists  a  sensitive  dependence  between  ply  pattern 
design  and  number  of  plies.  For  xample,  ply  patterns  to  con¬ 
struct  the  same  body  as  in  Figure  9  with  the  same  values  of  c 
and  T,  but  corresponding  to  N  =  90  and  91  are  shown  in  Figures 
10  and  11,  respectively.  This  relationship  is  also  accompanied 
by  marked  variations  in  tilt  angle  y.  While  minimal  variation 
occurs  within  a  given  body,  the  average  values  of  y  for  the 
cases  N  =  90,  91  are  approximately  17.6°,  15.2°,  respectively. 

The  equations  that  define  the  general  solution  for  the  EIS 
imply  certain  inequalities  that  must  be  satisfied  by  the  EIS 
parameters  to  guarantee  a  geometrically  admissible  solution. 

For  example,  because  of  Eq.  (23),  the  constant  c  must  be  equal 


to  or  less  than  the  smallest  value  of  r  in  the  body  of  revolu¬ 
tion  being  generated.  In  development  space,  the  angle  @  must  be 
single-valued,  so  that  Eqs.  (60)  and  (58)  require 


z _ -  z_  <  2ttc  /l  +  A2 

max  p - p, - 


(65) 


where  zmax  is  the  largest  value  of  z  in  the  body  of  revolution. 
Furthermore,  the  parameters  c,  T,  and  N  must  satisfy  Eq.  (41) . 
Finally,  fabrication  constraints,  which  can  limit  the  maximum 
and/or  minimum  size  of  a  body,  and  can  impose  restrictions  on 
central  angle  0,  must  also  be  taken  into  account  in  the  design 
of  involute  bodies. 

We  have  remarked  earlier  that  cylindrical  bodies  can  be 
constructed  via  non-cylindrical  involute  surfaces.  This  has  been 
illustrated  in  the  forming  of  the  cylindrical  region  of  the  body 
shown  in  Figure  9.  Conversely,  we  can  show  that  non-cylindrical 
bodies  can  be  formed  via  cylindrical  involute  surfaces.  This 
may  be  illustrated  by  invoking  Eq.  (49)  to  define  c  and  generating 
the  ply  pattern  corresponding  to  Figure  7  according  to  the 
present  model.  All  of  our  equations  remain  valid,  with  (27) , 

(32),  (35),  and  (38)  being  subject  to  the  interpretation  dis¬ 
cussed  in  connection  with  Eq.  (48).  For  values  of  the  parameters, 
N  =  200  and  T  =  .01,  the  ply  pattern  shown  in  Figure  12  is 
generated. 


DEBULKING  KINEMATICS 

In  practice,  involute  bodies  are  fabricated  from  sheets  of 
cloth-reinforced  composite  material.  In  its  initial  (as  received) 
form,  this  pre impregnated  cloth  material  is  uncured  and  relatively 


thick  (approximately  twice  the  final  cured  thickness) .  It  is 
in  this  form  that  the  material  is  cut  into  the  ply  pattern  and 
laminated  against  specially  designed  tooling.  Because  of  the 
large  thickness  in  this  form,  the  involute  surface  differs 
considerably  from  that  existing  in  the  final  state,  which  has 
governed  the  ply  pattern  design.  The  system  is  then  debulked 
under  a  slight  temperature  rise  and  the  application  of  pressuri¬ 
zation  on  either  the  inner  or  outer  surface,  with  the  opposite 
surface  being  controlled  by  tooling.  Subsequent  curing,  carbon¬ 
ization,  and/or  graphitization  establish  the  final  dimensions 
of  the  billet. 

Because  of  the  above  fabrication  procedure,  it  is  important 
to  define  the  configuration  of  the  involute  surface  corresponding 
to  material  in  its  initial,  undebulked  form.  It  is  also  necessary 
in  the  execution  of  the  debulking  phase  to  pressurize  against 
tooling  that  controls  the  profile  of  either  final  surface,  i.e., 

R^(z),  or  R 2 (z)  in  Figure  1.  It  is  also  desirable  that  the 
latter  ply  boundary  remains  in  contact  with  the  tooling  throughout 
the  process.  Unfortunately,  the  edge  of  the  ply  pattern  in  the 
initial  (undebulked)  configuration  cannot,  in  general,  coincide 
with  the  same  control  surface  as  that  in  the  final  configuration. 

An  exact  solution  for  the  initial  configuration  is  possible  if 

we  abandon  the  control  surface  requirement.  Thus,  we  simply 

formulate  the  inverse  problem:  given  the  ply  pattern,  N,  and  a 

new  ply  thickness,  define  the  corresponding  EIS.  In  this  case,  | 

we  assume  that  the  ply  pattern  construction  via  the  involute 

t« 
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method  (Figure  6)  remains  unchanged  while  a  new  value  of  c 


is  computed  by  use  of  Eqs.  (56)  and  (40),  i.e.. 


&'  ■  *  m 


where  primes  denote  undebulked  parameters.  Equations  (59)  and 
(62)  are  now  solved  with  the  aid  of  (56)  to  give  the  corresponding 
r,  2  coordinates. 


z  '  -  z  •  = 

P 


cot  X 


(r')2  =  c'R0cot2X  +  (c1)2  (68) 

while  a  single  value  of  z 1  (say  z£)  can  be  prescribed  arbitrarily. 

While  the  approach  outlined  in  the  previous  paragraph  is 
exact,  it  is  not  suitable  for  modern  involute  exit  cone  manufacture 


Therefore,  we  shall  now  consider  an  alternate  method  which  permits 
the  appropriate  edge  of  the  ply  peittern  to  lie  extremely  close  to 
the  control  surface  over  a  significant  region  in  the  initial  con¬ 
figuration  provided  the  ply  pattern  boundary  has  no  discontinuities 
in  the  region.  In  this  approach,  the  two  end  points  of  the  ply 
pattern  boundary  in  the  (smooth)  region,  which  we  shall  call 
reference  points,  lie  on  the  control  surface.  The  approach 
proceeds  as  follows : 


a.  Determine  the  ply  pattern  corresponding  to  the  final 
(cured)  state  in  the  usual  manner  and  compute  the  distance  L 
between  the  two  reference  points  in  the  R  plane.  We  let 

the  initial  position  of  the  first  point  coincide  with  its  final 
location,  while  the  second  point  is  merely  constrained  to  lie  on 
the  control  surface. 


b.  For  the  first  approximation,  let  the  initial  r,  6,  z 
coordinates  of  the  second  point  coincide  with  their  final  values. 
Using  Eqs.  (45)  and  (23),  compute  the  value  of  the  constant  c' 
for  the  initial  configuration. 

c.  Compute  the  initial  mapped  coordinates  (R* ,  @  ')  of  the 
first  reference  point  by  use  of  Eqs.  (56) -(59)  and  (62)  directly. 
Determine  the  initial  spatial  (r,  z)  position  of  the  second  point 
by  using  the  facts  that  it  must  lie  on  the  control  surface  and 
its  image  is  at  distance  L  from  that  of  the  first  point  in  the 
ply  pattern.  Now  compute  the  (R1,  ©  ')  coordinates  of  the 
second  point.  It  should  be  noted  that  the  '  coordinate 
involves  an  arbitrary  (non-essential)  constant  z^. 

d.  At  this  time,  we  have  the  mapped  coordinates  of  each 
reference  point  in  two  coordinate  systems,  (R,  (h)  )  and  (R',  (h)  ' ) 
Since  the  two  images  of  each  reference  point  must  coincide,  we 
now  establish  the  relationship  between  the  two  sets  of  coordinates 

e.  By  using  the  previously  derived  relationship,  determine 
the  (R' ,  (S)  ' )  coordinates  of  the  ply  pattern.  Finally,  apply 
the  inverted  form  of  Eqs.  (56)  —  (59)  and  (62),  along  with  (43) 
and  (23) ,  to  define  the  spatial  positions  of  points  on  the  ply 
pattern  boundary. 

Adjustment  of  the  constant  c'  can  be  accomplished  to  achieve 
a  closer  fit  between  the  initial  ply  pattern  boundary  and  the 
control  surface,  whereupon  steps  (a)  through  (e)  are  repeated. 

The  above  procedure  involves  only  elementary  geometric  relations 
in  addition  to  repetitious  use  of  Eqs.  (56)-(59),  (62),  (45),  and 
(23).  Hence,  these  equations  are  not  shown  here.  Use  of  this 
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procedure  (first  approximation  only)  leads  to  the  profile  shown 
in  Figure  13  for  the  ply  pattern  of  Figure  8,  where  the  profile 
is  based  upon  a  ply  thickness  of  .02  in,  while  the  reference 
points  referred  to  in  the  approach  are  selected  as  the  extremities 
of  the  inner  conical  surface.  All  points  on  the  ply  boundary  do 
not  contact  this  surface  although  the  deviations  are  imperceptible 
at  the  scale  of  Figure  13. 

CONCLUDING  REMARKS 

The  general  problem  of  involute  construction  -  formation  of 
a  body  of  revolution  by  laminating  identical  sheets  of  constant 
thickness  -  has  been  treated  in  detail.  The  two  transformations 
that  were  applied  -  movement  of  each  point  of  the  solution 
surface  a  constant  distance  in  the  surface  normal  direction 
followed  by  a  rigid  body  rotation  of  the  surface  -  guarantee 
complete  filling  of  the  required  volume  with  no  gaps  or  overlapping 
regions.  The  exact  involute  surface  (EIS)  thus  satisfier  all 
requirements  provided  the  profile  of  the  body  of  revolution  is 
piecewise-continuous.  The  question  of  the  uniqueness  of  the 
particular  pair  of  transformations  that  were  applied  to  carry  the 
EIS  into  itself  has  not  been  addressed.  However,  the  solution 
that  satisfies  the  particular  transformations  and  input  data  (47).' 
i.e.,  Eqs.  (45)  and  (23),  does  represent  a  unique  surface,  and, 
of  course,  satisfies  the  problem  requirements.  Thus,  the  very 
important  question  of  existence  of  a  solution  to  the  general 
problem  of  involute  construction  has  been  settled. 

We  have  rigorously  shown  that  the  EIS  assumes  the  geometric 
form  of  a  helical  convolute,  which  was  evidently  first  postulated 


in  unpublished  work  by  Savage  (Ref.  7) .  The  helical  convolute 
is  a  surface  of  single  curvature  that  is  capable  of  development 
by  standard  methods  of  descriptive  geometry  (Ref.  8) .  The  closed- 
form  equations  that  describe  this  developed  surface,  called  the 
ply  pattern,  have  also  been  presented  here.  These  are  Eqs. 

( 56 ) — ( 59 ) ,  (62),  and  (64).  Certain  inequalities  that  must  be 
satisfied  by  the  EIS  parameters,  including  Eqs.  (41)  and  (65) , 
have  also  been  established. 

The  existence  of  a  developable  solution  surface  implies  that 
a  defect-free  body  is  a  potential  product  of  this  work.  This  is 
in  sharp  contrast  to  the  severe  fabric  distortion  and  wrinkles 
observed  (Ref.  9)  in  contemporary  involute  exit  cones.  Demon¬ 
stration  of  this  point  and  some  validation  of  the  mathematical 
formulation  have  been  provided  by  the  construction  of  the 
experimental  model  shown  in  Figure  9 . 

Finally,  an  approach  to  describe  the  kinematics  of  the 
debulking  process  has  been  established.  By  use  of  this  model, 
one  can  approximate  quite  accurately  the  condition  that  a  ply 
boundary  remains  in  contact  with  a  fabrication  tooling  surface, 
provided  the  ply  boundary  contains  no  discontinuity  in  the 
region.  Another  use  of  this  model  is  to  predict  the  displace¬ 
ment  field  as  a  function  of  instantaneous  ply  thickness  in  the 
debulking  process. 


APPENDIX  A 

VOLUME  AND  SURFACE  AREA 


Letting  S  represent  the  area  of  the  EIS  and  referring 
to  Figure  4,  we  observe  that 


dS  =  ds  — -  sin  (A' 

cos  y 

where  ds  is  the  infinitesimal  arc  length  of  an  involute  in  the 
plane  z  =  constant.  From  Reference  4,  ds  is  given  by 


ds 


-dr 
sin  a 


(A 


Hence,  (A-l)  becomes 


ds  =  ~sin  \p  dr  dz 
cos  y  sin  a 


(A- 


Substitution  of  Eqs.  (7) ,  (23) ,  and  (46)  and  use  of  elementary 
trigonometry  brings  (A- 3)  to 


ds  =  ~cd  +  A*)  **  cosot  dadz 
.  sin3  a 


Integrating  and  substituting  Eqs.  (23)  and  (40)  yields 

S  =  ST  J  lR!(z)  "  Ri(z)]  dz 


(A- 


(A- 


where  the  nomenclature  of  Figure  1  has  been  used.  Hence,  the 
known  result 

V  =  NTS  (A' 

is  identically  satisfied  by  the  present  formulation.  Here,  V 
stands  for  the  volume  of  the  body  of  revolution. 
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APPENDIX  B 
HELICAL  ANGLE 


An  important  parameter  in  the  analysis  of  mechanical  response 
of  involute  exit  cones  is  the  helical  angle  <j>  (Ref.  6)  ,  defined 
as  the  local  angle  between  the  warp  fibers  in  a  ply  and  the 
meridian.  Referring  to  Figure  14,  we  have  the  relation 

♦  =  <DW  -  ®t  <b-i> 

where  ®w  is  the  constant  angle  between  the  warp  direction  in 
the  ply  pattern  and  an  arbitrary  axis  and  ®t  gives  the  orienta¬ 
tion  of  the  mapped  meridian  (denoted  by  t  in  Figure  14) . 

In  Figure  14,  the  dashed  line  marked  s  represents  the  mapped 
position  of  the  local  vector  s.  Observe  that  the  angle  i|>,  given 
by  Eq.  (7) ,  is  preserved  in  the  mapping.  Hence,  we  see  that 

®t  =  ir/2  -  ij>  +  X  +  (h)  (B-2) 

Hence  by  use  of  Eqs.  (57)  and  (61)-(63),  we  arrive  at 

(Dv.  =  7T  -  i(i  +  - - -  [§(z  -  z_)  -  cota]  (B-3) 

VT777  p 

which,  along  with  (B-l) ,  produces  the  desired  result. 
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APPENDIX  C 

FIGURES 


2.  Lxact  Involute  Surface  of  a 
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